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Abstract — In this paper, fast-decodable lattice code construc- 
tions are designed for tlie nonorthogonal amplify-and-forward 
(NAF) multiple-input multiple-output (MIMO) channel. The con- 
structions are based on different types of algebraic structures, e.g. 
quaternion division algebras. When satisfying certain properties, 
these algebras provide us with codes whose structure naturally 
reduces the decoding complexity. The complexity can be further 
reduced by shortening the block length, i.e., by considering 
rectangular codes called less than minimum delay (LMD) codes. 

I. Introduction 

The quaUty of wireless long distance communications can 
be significantly improved by using cooperative diversity tech- 
niques. Cooperating relays can be positioned between the 
source station and the destination to aid the transmission 
by either amplifying and forwarding (AF) or decoding and 
forwarding (DF) the signal. Spatially separated terminals will 
allow an increment in the diversity in a distributed manner. 
Depending on the application, a one-hop or multi-hop trans- 
mission is called for. Here, we consider multi-hop distributed 
space-time codes employing a half-duplex NAF protocol H], 
121. It is known 111 that the NAF protocol outperforms all other 
AF protocols since, as opposed to orthogonal protocols, it can 
keep transmitting also during the transmission of the relays. 
In addition, the AF protocols are less complex than the DF 
protocols. This type of low cost relay systems are called for 
in e.g. digital video broadcasting (DVB) |4J. 

In Is) and ||5|, Yang et al. and Hollanti et al. proposed block- 
diagonal space-time code constructions for the asymmetric 
MIMO channel with or without relays. The constructions arise 
from cycUc division algebras constructed over a higher degree 
center A nonvanishing determinant (NVD) is then achieved 
by forming a block-diagonal matrix consisting of the left 
regular representation of the algebra and its Galois conjugates 
from the center to the base field. It was also shown 13] that 
a block-diagonal structure together with the NVD property 
is enough to achieve the diversity-multiplexing gain tradeoff 
(DMT) also in the asymmetric case, where the number of 
transmit antennas is strictly bigger than the number of receive 
antennas, and hence the corresponding lattice is not full. 
Motivated by this and the urge for complexity reduction of 
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MIMO codes in general, we impose further properties that 
the algebras and the constructions should satisfy in order to 
reduce the complexity. Our study reveals a trade-off between 
the coding gain and decoding complexity. Related work has 
been carried out by, among others, Rajan et al. (see e.g. |fT2l). 
They considered fast-decodable distributed space-time codes 
arising from Clifford algebras. Our work differs from theirs 
in that our codes achieve the NVD property and hence the 
asymmetric DMT. The codes proposed in this paper moreover 
have a nice algebraic structure which makes analyzing the 
codes easier 

List of contributions: 

• Explicit fast-decodable space-time relay codes are pro- 
posed. 

• All the codes have full diversity, some even NVD. To the 
best of the authors' knowledge, these are the first fast- 
decodable distributed space-time codes with NVD. 

• The constructions arise from quaternion or other type of 
algebras making it easy to determine the coding gain, 
complexity and other properties of the codes. 

• Although our explicit examples are for the case when the 
source and the relays each have only one antenna, the 
constructions are easily generalizable to other numbers 
of antennas and relays as well. 

• We demonstrate a performance-complexity tradeoff re- 
sulting from the used method. 

• Finally, we analyze the worst-case decoding complexity 
of the proposed codes. 

Let us finish this introductory section by giving a couple 
essential definitions. 

Definition 1.1: If the code C consisting of matrices X 
satisfies 

min ActlX^X) > k > 0, 

we say that C has the nonvanishing determinant property 
(NVD). 

In case of square matrices, we simply refer to dci{X) when 
talking about NVD. 

There are multiple definitions of rate, but we will consis- 
tently use the following. 



Definition 1.2: Let Bi, . . . ,Bk G Af„jxT(C) be the gener- 
ator matrices (over R) of a rank k code C, so 



C = ^B,g„ 



where e Z, e.g. PAM symbols. The rate R of the code is 
then 

R = fc/T 

(real) dimensions per channel use (dpcu). 

Note that the commonly used rate in complex dimensions 
per channel use is R/2 when using the above notation. 

II. System model for the NAF relay channel 

For ease of notation, we only define the single-relay model, 
the generalization to multi-hop is straightforward. Following 
||3|, let us denote by Xi the signals transmitted from the source, 
and by Yr the signal received by the relay which the relay then 
amplifies and forwards as BY,-. The number of relays and 
the number of antennas at the source, relays and destination 
are denoted by N,ns,nr,nd, respectively. We assume rir is 
the same for all relays r = 1,...,A^. To be realistic, we 
assume Us > rir- The destination is observing Yi and Y2 in 
consecutive time instances, and we have 



Yi 

Y2 



v/tti SNRFXi + Vi 

^/iTipSNRHXi + Wi 

yJ-K-i SNRG{BYr) + ^f^^SNRFX2 + V2. 



where Vi, W are the additive white gaussian noise matrices 
and H, G are the Rayleigh distributed channel matrices. 
The power allocation tt^ factors are chosen so that SNR 
denotes the received SNR per receive antenna at the destina- 
tion. We assume perfect channel state information (CSI) at the 
receivers, while the transmitters have none. For more details, 
we refer to fS). 

III. On fast decodability 

Maximum-likelihood decoding amounts to searching the 
code C for the codeword 



Z = arg'min{\\Y - HX\\\}xeC-, 



(1) 



closest to the received matrix Y with respect to the squared 
Frobenius norm. 

Consider a code C of Q-rank k, i.e., each codeword X 
is a linear combination X^iLi Bigi of generating matrices 
_Bi, . . . , Bk, weighted by coefficients gi, . . . ,gk, which are 
PAM information symbols. The matrices Bi, . . . ,Bk therefore 
define our code. Each n,- x T matrix HBi corresponds, via 
vectorization, to a vector b,; e R^^"'' obtained by stacking the 
columns followed by separating the real and imaginary parts 
of HBi. We define the (generating) matrix 

S= (bi,b2,...,bfc) e M2TrvxfcW, 



SO every received codeword can be represented as a real vector 
Bg, with g = [gi, . . . ,gk)^ having coefficients in the real 
alphabet S in use. 

Now finding argmin{\\Y — HX\\y}xec becomes equiv- 
alent to finding argmin{\\y — -Bg|||;}gg|5|fc with respect to 
Euclidean norm, where y is the vectorization of the received 
matrix Y . The latter search is performed using a real sphere 
decoder ||6l, with the complexity of exhaustive search amount- 
ing to \S\^, as the coefficients of g run over all the values 
of S. The complexity of decoding can, however, be reduced 
if the code has additional structure ff\. Performing a QR 
decomposition of B, B = QR, with Q^Q — I, reduces finding 
argminlWy — i3g|||;}g to minimizing 



\y-QRs\\l^\\Q^y-Rs\\l, 



(2) 



where R is an upper right triangular matrix. The number and 
structure of zeros of the matrix R may improve the decoding 
complexity (formally defined HI to be the minimum number of 
vectors g over which the difference in (|2]l must be computed). 
When the structure of the code allows for the degree (i.e., the 
exponent of |5|) of decoding complexity to be less than the 
rank of the code, we say that the code is fast-decodable. 

More precisely, we have the following definitions from Q. 

Definition 3.1: A space-time code is said to be fast- 
decodable if its R matrix has the following form: 



R = 



A 




Bi 

i?2 



where A is a diagonal matrix and R2 is upper-triangular 

The authors of |7i give criteria when the zero structure of 
R coincides with that of M, where M is a matrix capturing 
information about orthogonality relations of the basis elements 
of Bf. 



Ml 



k,l 



\BlBi+BjBu\ 



F- 



(3) 



In particular, 1 3 Lemma 2] shows that if M has the structure 



M 



A 

B2 



Bi 
B3 



A 




Bi 
Ri 



where A is diagonal, then R 

We could thus rephrase Defnition 13. II in terms of M. 

Next we recall the class of codes which allows groups of 
symbols to be decoded independently of one another. 

Definition 3.2: A space-time code of dimension K is called 
g-group decodable if there exists a partition of {1, ... , A'} into 
g nonempty subsets jTi, . . . , Jg, so that the matrix A/(,fc = 
when I, k are in disjoint subsets Ji, Jj. 

In this case, as shown in Q, the matrix R has the form 
~ 
R= ••• 



where each Ri is a square upper triangular matrix. Hence, 
the symbols Xk and xi can be decoded independently when 
their corresponding basis matrices Bk and Bi belong to 
disjoint subsets of the partition. 

Remark 1: Note that a simple computation shows that the 
zero structure of M is stable under premultiplication of Bi by 
a channel matrix H. In general, the same does not hold for 



R. 

By the above discussion, in order to demonstrate fast- 
decodability (resp. g-group decodability), it suffices to find 
an ordering on the basis elements Bi, which results in the 
desired zero structure of M. We proceed to do that for the 
proposed relay codes. 



IV. Minimum delay codes for N 

RELAYS 



2 AND N 



We demonstrate that the codes we obtain are conditionally 
4-group decodable. Recall from ITTl that a code is called 
conditionally g-group decodable if there exists a partition of 
the indices {!,... ,K} of basis elements into g + I disjoint 
subsets J'l, .. ., Jn, J"^ such that 



\BlB,r, + bIBiWf = yie J^ym e Jj,i ^ j. 



In this case, the sphere decoding complexity order reduces 



to IS*! 



|J"-|+maxi<i<j, \Ji\ 



A. Virtual 6x2 MIMO channel with N = 3,ns = rir = 
1, rid = 2 



We proceed to show the rate four (4 PAM symbols per 
channel use) relay construction consisting of 6 x 6 matrices. 



Proposition 4.1: Define the code 

[X 

C - {ar{X)} 





t{X) 
T^X) 



where X is a matrix of the form 

^ = {vTid 

with c,dE Z{i, (i), a : —i. 



lla{d) 
a{c) 



Then the code C is of rank 24 and (real) decoding complex- 
ity 151^^, and has the NVD property. 



Proof: Define K = QiV^), K' 
Ct, then K' = K{C + C-^)■ 



^). Let ( denote 



(-11,-1) 



K' 




Q(C7) = K' 



First we note that the algebra (—11,-1 



is division. 



This follows from techniques of lfT3l : we apply ifTSl Theorem 
7.1] while noting that Fna contains no element of order 4, 
i.e.,— 1 is not a square in Fus, which is the residue field of 
the prime 11 in ^(Ct)- 

Now note that (-11,-1)a" = (-1,-11)^', so let Q = 
(—1, — 11)/^/. After conjugation that does not affect the deter- 
minant but does aid energy balance and decoding complexity, 
it has the following A''-basis: 







'1 




'i ■ 








1 




-i 









s/Tli 







-^/T^ 


93 = 







,94 = 








This means that Q is generated over K 
following 12 matrices 



7) by the 



for i = 1, ... ,4. Extend this to a Q-basis by letting r^.j = 
^/—^Ti-4^ j for i = 5, . . . , 8. Then a Z-basis of C can be given 
by 



{0tT{Yl.j)}l<d. 



<8,j<3 



(4) 



and is of size 24. Indeed, the rank of C is 24, since each code- 
word X is an element of quaternion algebra (—1,-11)^^', and 
hence encodes 4 symbols from K' = Q(C7), or equivalently 
24 real symbols. 

Now let T : (^7 H-> be a generator of Gal{K' /K). When 
the coefficients of codewords are algebraic integers, the code is 
NVD. This follows from the fact that the determinant of each 
codeword is fixed by both r and a, hence it is an element 
of (Q(a/^7) = K. Moreover it is nonzero whenever X is 



nonzero, since Q was shown to be division. 

We show that C is conditionally 4-group decodable with 
complexity \S\^^^\ conditioned on decoding symbols corre- 
sponding to {r.s i, . . . , Fg^a}, the complexity of decoding 
symbols corresponding to {Fi i, . . . , r4_3} is at most |5'|'^, 
where S is the underlying alphabet. For that, note that when 
A = Fj j-, B = Fi' jv, for all j, j' and for i ^ i' , we have 

AB'^ + BA'^ = 0. 

Same follows for ar{A),ar{B), i.e., we have: 

ar{A)ar{B)'^ + = 0. 

Let F = [Q!r(Fi,i), . . . ,a^(F8,3)] be the Ust of 24 genera- 
tors of C from (|4|i in lexicographical order Then the matrix 
AI = Alij from Equation (O capturing orthogonality relations 
on F has the following structure: 



M 



* 

* 












(5) 



where each coefficient of the matrix above is a 3 x 3 matrix, 
which is when the coefficient is 0. 

■ 

B. Virtual 4x2 MIMO channel with N = 2,ns ~ 71,. = 1, 
rid = 2 

We use a similar idea as in the 6x6 case to construct fast- 
decodable rate four relay codes consisting of 4 x 4 matrices. 

Proposition 4.2: Let the code C consist of matrices 

^ _ c -V2(T(d) 

t{X) ^ ^ V2d aic) 



X 




with c,d G Ok, a : y/a i— > — Again we have conjugated 
the original matrix in order to aid decodabiUty and energy 
balance. The minimum determinant is invariant under such 
conjugation. 

Each matrix X represents an element from the quaternion 
algebra (0,7);^ over a biquadratic field K = Q{i, ,Jp). Then 
for values a = 5,7 = — 2,p = 31 the resulting code is a 
fully diverse NVD code of rank 16. It is conditionally 4-group 
decodable with (real) decoding complexity \S\^^. 

Proof: We use similar techniques to the previous proof 
to establish that (a, is division. First we establish that the 
prime ideal 5Z splits completely in K. Using the fact that its 
residue field is isomorphic to F5, we conclude that 7 is not 
a square modulo 5 in the integers of K. Hence using ifTsl 
Theorem 7.1], we conclude that Q = {a,^)K is division. The 
generators of Q over Q(i) are 



We use the image of these generators under ar to generate 
the code C over Z. To check fast decodability, we verify the 
relations 



ar{X)^ar{Y) + ar{Y)'^ar{X) 



0, 



X = F,j, r = Fj. J/, where 1 < i 7^ i' < 4 and 1 < j, j' < 3. 

The ordering ^^(Fia), . . . ,a^(iFi4), . . . ,a7-(iF4_2) gives 
the matrix M with the same zero structure as in Equation Q 
only now each coefficient is a 2 x 2 matrix. Hence the code is 
conditionally 4-group decodable with complexity 15*1^°. More 
precisely, conditioned on decoding symbols corresponding to 
{aT-(iFi_i), . . . , ai-(iF4 2)}, the complexity of decoding the 
rest of the symbols is at most \S\'^. 

■ 

Remark 2: Here we have concentrated on a real sphere de- 
coding process. Note, however, that Fi.i, F1.2, • . • , F4 1, F4 2 
gives a Z(i)-basis of C. We can verify the relations 

ar{X)'far{Y) + ar{Y)'far{X) = 0, 

X ^ Fj.j , r = Fj/j/, where 1 < i ^ i' < A. Hence the 
obtained code is 4-group decodable, of complexity degree 2 
when using a complex decoder. 

V. Less than minimum delay codes for N = 2 

RELAYS 

In order to further reduce the complexity, we will shorten 
the block length. Such less than minimum delay codes have 
been considered in |[9l • To this end, let us start by constructing 
a code with dimension rate R = 2, i.e., the lattice is of rank 
4 and the code matrix transmits two real dimensions (e.g. two 
PAM symbols) per channel use. 

A. Virtual 4x1 MIMO channel with N = 2,ns = rir = 
1, nd = l 

Let us consider the 8th cyclotomic extension (Q(C8)/Q, and 
denote its Galois group by {1, t : \/2 t-^ — \/2,* , t*}, where 
* denotes complex conjugation. Our code will simply consist 
of matrices 



X{ai,a2) 



( 



X 

X* 









t(x) 
t(x 



\ 



F,; 



1, 



,4. 



where x ~ ai + a2C8 and o-i G ^[*]- Due to the fact that the 
code matrix only contains two QAM symbols, the complexity 
will automatically be at most \S\'^, where S is the underlying 
real alphabet. However, if we do a smart ordering of the basis 
elements as {Bi = X(l, 0), X{i, 0), X(0, 1), X(0, i)}, the 4x 
4 matrix R will have the form described in Definition 13 . 1 1 with 
a 2 X 2 matrix A. This is due to the fact that whenever we 
have a totally real basis element Bi and a totally imaginary 
one Sj, this will result in mT[{HBi)^ HBj] ^ 0. Thus, we 
have reduced the worst case complexity to |S'p, that is, by 
62.5% compared to the complexity \S\^ of a general square 
code with the same rate. 



Proposition 5.1: The matrices X'^X, X ^ 0, have the 
NVD property. 

Proof: The proof is straightforward. Namely, 
det(Xl'X) = \xt\x)x* t(x)* = i^(^^^)q,(x) G 4Z, 
and hence dct(Xl'X) > 4. ■ 

Next, we extend the above construction to the case = 2 
ideally calling for a rate four code. 

B. Virtual Ax 2 MIMO channel with N ~ 2,ns = rir = 
1, nd = 2 

Let us next construct a rank 8 lattice in order to have higher 
multiplexing of 4 real dimensions per channel use. We start 
by adjoining \/5 to the above extension, i.e., we consider 
QCCs; V5)/Q and denote the corresponding maps fixing Q(i) 
by {1, T : V2 i-^ -V2, r : ^/5 ^ - ^5, rr}. 

The code matrix now looks like 

r{iyx) 
A(ai,a2,a3,a4j = g ^^^^^ 

y Tr{vx) J 

where x = oi + a2C8 + + 0,4(^6, 9 ~ '^^■/^ , cLi G "^[i], 
and J/ — \ + i — iO generates a principal ideal that will 
make the code lattice orthogonal. This field extension is the 
same as the one used for the extended golden algebra in m. 
The complexity of the code is at most \S\^ . Similarly to the 
Ud = 1 case, this can be further reduced to by ordering 
the basis as {Bi = 0, 0, 0), B2 = X{i, 0, 0, 0), . . .}. We 
have reduced the complexity by 56.25% compared to a general 
square code with the same rate and complexity \S\^^. 

Proposition 5.2: The matrices X^X, X ^ 0, have full 
rank, i.e., a code consisting of the matrices X has full diversity. 

Proof: Again, the proof is very simple. We have 
det(XtX) = + \r{iyx)\^){\T{iyx)\^ + \Tr{iyx)\^) > 0. 

■ 

Remark 3: We want to point out that, unfortunately, there 
is no free lunch. Namely with the above construction method 
increasing the code rate causes a degradation in the coding 
gain. At the same time, reducing the delay from four to 
two channel uses will naturally lower the maximum rank 
(and diversity) the code matrix can achieve when compared 
to a square code matrix. Hence, we observe a performance- 
complexity tradeoff implying that while we can indeed reduce 
the decoding complexity by reducing the code length, we are 
likely to face slightly worse performance caused by the reduc- 
tion in diversity. Increasing the rate (performance) will also 
here require dropping the NVD requirement, again indicating 
a tradeoff type behavior 

VI. Conclusions and future work 

We proposed explicit relay codes with a fast-decodable 
structure and NVD for different number of antennas 
and relays. It was shown that the method used implies 
a performance-complexity tradeoff. In other words, fast- 
decodable codes with NVD were proposed, while at the same 
time it was noted that relaxing on the NVD property allows 



for further complexity reductions. One efficient way to reduce 
the complexity is to employ less than minimum delay codes 
that by construction already halve the complexity compared 
to general minimum delay codes. It remains to be investigated 
how the codes perform compared to other distributed codes 
that either have higher complexity or lack NVD. Related 
results lfT4l will be posted to arXiv in near future. 
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